Multi-parametric MRI-based machine learning model for prediction of WHO grading in patients with meningiomas.

OBJECTIVE: The purpose of this study was to develop and validate a nomogram combined multiparametric MRI and clinical indicators for identifying the WHO grade of meningioma. MATERIALS AND METHODS: Five hundred and sixty-eight patients were included in this study, who were diagnosed pathologically as having meningiomas. Firstly, radiomics features were extracted from CE-T1, T2, and 1-cm-thick tumor-to-brain interface (BTI) images. Then, difference analysis and the least absolute shrinkage and selection operator were orderly used to select the most representative features. Next, the support vector machine algorithm was conducted to predict the WHO grade of meningioma. Furthermore, a nomogram incorporated radiomics features and valuable clinical indicators was constructed by logistic regression. The performance of the nomogram was assessed by calibration and clinical effectiveness, as well as internal validation. RESULTS: Peritumoral edema volume and gender are independent risk factors for predicting meningioma grade. The multiparametric MRI features incorporating CE-T1, T2, and BTI features showed the higher performance for prediction of meningioma grade with a pooled AUC = 0.885 (95% CI, 0.821-0.946) and 0.860 (95% CI, 0.788-0.923) in the training and test groups, respectively. Then, a nomogram with a pooled AUC = 0.912 (95% CI, 0.876-0.961), combined radiomics score, peritumoral edema volume, and gender improved diagnostic performance compared to radiomics model or clinical model and showed good calibration as the true results. Moreover, decision curve analysis demonstrated satisfactory clinical effectiveness of the proposed nomogram. CONCLUSIONS: A novel nomogram is simple yet effective in differentiating WHO grades of meningioma and thus can be used in patients with meningiomas. CLINICAL RELEVANCE STATEMENT: We proposed a nomogram that included clinical indicators and multi-parameter radiomics features, which can accurately, objectively, and non-invasively differentiate WHO grading of meningioma and thus can be used in clinical work. KEY POINTS: • The study combined radiomics features and clinical indicators for objectively predicting the meningioma grade. • The model with CE-T1 + T2 + brain-to-tumor interface features demonstrated the best predictive performance by investigating seven different radiomics models. • The nomogram potentially has clinical applications in distinguishing high-grade and low-grade meningiomas.

Comprehensive analysis of cuproptosis-related long non-coding RNAs in prognosis, immune microenvironment infiltration and chemotherapy response of hepatocellular carcinoma.

The objective of this study is to explore the relationship between cuproptosis-related long noncoding RNAs (lncRNAs) in hepatocellular carcinoma (HCC). RNA-seq data, including lncRNAs and related clinical information of HCC patients, were downloaded from The Cancer Genome Atlas database. A signature composed 3 cuproptosis-related lncRNAs was constructed by LASSO analysis, and HCC patients were classified into high- and low-risk groups. Patients in the high-risk group had a poorer prognosis compared with the low-risk group. Univariate Cox and multivariate Cox regression analyses confirmed that the signature model was an independent risk factor compared to other clinical biomarkers. Furthermore, gene set enrichment analysis indicated that metabolism-related pathways were enriched in low-risk group, including drug metabolism, and fatty acid metabolism. Further research demonstrated that there were markedly differences in drug response between the high- and low-risk group. Immune related analysis showed that the most type of immune cells and immunological function in the high-risk group were different with the risk-group. Finally, TP53 mutation rate and the tumor mutational burden in the high-risk group were higher compared with the low-risk group. In conclusion, we constructed a prognostic signature based on the expression of cuproptosis-related lncRNAs to predict HCC patients' prognosis, drug response and immune microenvironment, and further research will be conducted to uncover the mechanisms.

Integrative Analysis of Inflammatory Response-Related Gene for Predicting Prognosis and Immunotherapy in Glioma.

Inflammatory response plays a crucial role in the development and progression of gliomas. Whereas the prognostic esteem of inflammatory response-related genes has never been comprehensively explored in glioma, the RNA-seq information and clinical data of patients with glioma were extracted from TCGA, CGGA, and Rembrandt databases. The differentially expressed genes (DEGs) were picked out between glioma tissue and non-tumor brain tissue (NBT). Then, the least absolute shrinkage and selection operator (LASSO) regression analysis was performed to construct the prognostic signature in the TCGA cohort and verified in other cohorts. Kaplan-Meier survival analyses were conducted to compare the overall survival (OS) between the high and low-risk groups. Univariate and multivariate Cox analyses were subsequently used to confirm the independent prognostic factors of OS, and then, the nomogram was established based them. Furthermore, immune infiltration, immune checkpoints, and immunotherapy were also probed and compared between high and low-risk groups. The four genes were also analyzed by qRT-PCR, immunohistochemistry, and western blot trials between glioma tissue and NBT. The 39 DEGs were identified between glioma tissue and NBT, of which 31 genes are associated to the prognosis of glioma. The 8 optimal inflammatory response-related genes were selected to construct the prognostic inflammatory response-related signature (IRRS) through the LASSO regression. The effectiveness of the IRRS was verified in the TCGA, CGGA, and Rembrandt cohorts. Meanwhile, a nomogram with better accuracy was established to predict OS based on the independent prognostic factors. The IRRS was highly correlated with clinicopathological features, immune infiltration, and genomic alterations in glioma patients. In addition, four selective genes also verified the difference between glioma tissue and NBT. A novel prognostic signature was associated with the prognosis, immune infiltration, and immunotherapy effect in patients with gliomas. Thus, this study could provide a perspective for glioma prognosis and treatment.

Epigenetically silenced lncRNA SNAI3-AS1 promotes ferroptosis in glioma via perturbing the m6A-dependent recognition of Nrf2 mRNA mediated by SND1.

BACKGROUND: Ferroptosis has been linked to tumor progression and resistance to antineoplastic therapy. Long noncoding RNA (lncRNA) exerts a regulatory role in various biological processes of tumor cells, while the function and molecular mechanism of lncRNA in ferroptosis are yet to be clarified in glioma. METHODS: Both gain-of-function and loss-of-function experiments were employed to investigate the effects of SNAI3-AS1 on the tumorigenesis and ferroptosis susceptibility of glioma in vitro and in vivo. Bioinformatics analysis, Bisulfite sequencing PCR, RNA pull-down, RIP, MeRIP and dual-luciferase reporter assay were performed to explore the low expression mechanism of SNAI3-AS1 and the downstream mechanism of SNAI3-AS1 in ferroptosis susceptibility of glioma. RESULTS: We found that ferroptosis inducer erastin downregulates SNAI3-AS1 expression in glioma by increasing the DNA methylation level of SNAI3-AS1 promoter. SNAI3-AS1 functions as a tumor suppressor in glioma. Importantly, SNAI3-AS1 enhances the anti-tumor activity of erastin by promoting ferroptosis both in vitro and in vivo. Mechanistically, SNAI3-AS1 competitively binds to SND1 and perturbs the m6A-dependent recognition of Nrf2 mRNA 3'UTR by SND1, thereby reducing the mRNA stability of Nrf2. Rescue experiments confirmed that SND1 overexpression and silence can rescue the gain- and loss-of-function ferroptotic phenotypes of SNAI3-AS1, respectively. CONCLUSIONS: Our findings elucidate the effect and detailed mechanism of SNAI3-AS1/SND1/Nrf2 signalling axis in ferroptosis, and provide a theoretical support for inducing ferroptosis to improve glioma treatment.

TRIM56 acts through the IQGAP1-CDC42 signaling axis to promote glioma cell migration and invasion.

Diffuse invasion is an important factor leading to treatment resistance and a poor prognosis in gliomas. Herein, we found that expression of the tripartite motif containing 56 (TRIM56), a RING-finger domain containing E3 ubiquitin ligase, was markedly higher in glioma than in normal brain tissue, and was significantly correlated with malignant phenotypes and a poor prognosis. In vitro and in vivo experimental studies revealed that TRIM56 promoted the migration and invasion of glioma cells. Mechanistically, TRIM56 was transcriptionally regulated by SP1 and promoted the K48-K63-linked poly-ubiquitination transition of IQGAP1 at Lys-1230 by interacting with it, which in turn promoted CDC42 activation. This mechanism was confirmed to mediate glioma migration and invasion. In conclusion, our study provides insights into the mechanisms through which TRIM56 promotes glioma motility, i.e., by regulating IQGAP1 ubiquitination to promote CDC42 activation, which might be clinically targeted for the treatment of glioma.

A prognostic NAD+ metabolism-related gene signature for predicting response to immune checkpoint inhibitor in glioma.

Background: Nicotinamide adenine dinucleotide (NAD+) metabolism is involved in a series of cancer pathogenesis processes, and is considered a promising therapeutic target for cancer treatment. However, a comprehensive analysis of NAD+ metabolism events on immune regulation and cancer survival has not yet been conducted. Here, we constructed a prognostic NAD+ metabolism-related gene signature (NMRGS) associated with immune checkpoint inhibitor (ICI) efficacy in glioma. Methods: 40 NAD+ metabolism-related genes (NMRGs) were obtained from the Reactome database and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Glioma cases with transcriptome data and clinical information were obtained from Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA). NMRGS was constructed based on the calculated risk score using univariate analysis, Kaplan-Meier analysis, multivariate Cox regression, and nomogram. This NMRGS was verified in training (CGGA693) and validation (TCGA and CGGA325) cohorts. The immune characteristics, mutation profile, and response to ICI therapy were subsequently analyzed for different NMRGS subgroups. Results: Six NAD+ metabolism-related genes, including CD38, nicotinamide adenine dinucleotide kinase (NADK), nicotinate phosphoribosyltransferase (NAPRT), nicotinamide/nicotinic acid mononucleotide adenylyltransferase 3 (NMNAT3), poly(ADP-Ribose) polymerase family member 6 (PARP6), and poly(ADP-Ribose) polymerase family member 9 (PARP9), were ultimately used to construct a comprehensive risk model for glioma patients. Patients in the NMRGS-high group showed a poorer survival outcome than those in the NMRGS-low group. The area under curve (AUC) indicated that NMRGS has good potential in glioma prognostic prediction. A nomogram with improved accuracy was established based on independent prognostic factors (NMRGS score, 1p19q codeletion status, and WHO grade). Furthermore, patients in the NMRGS-high group showed a more immunosuppressive microenvironment, higher tumor mutation burden (TMB), higher human leucocyte antigen (HLA) expression and a more therapeutic response to ICI therapy. Conclusions: This study constructed a prognostic NAD+ metabolism-related signature associated with the immune landscape in glioma, which can be used for guiding individualized ICI therapy.

Exploring the multidimensional heterogeneities of glioblastoma multiforme based on sample-specific edge perturbation in gene interaction network.

Glioblastoma multiforme (GBM) is the most malignant brain cancer with great heterogeneities in many aspects, such as prognosis, clinicopathological features, immune landscapes, and immunotherapeutic responses. Considering that gene interaction network is relatively stable in a healthy state but widely perturbed in cancers, we sought to explore the multidimensional heterogeneities of GBM through evaluating the degree of network perturbations. The gene interaction network perturbations of GBM samples (TCGA cohort) and normal samples (GTEx database) were characterized by edge perturbations, which were quantized through evaluating the change in relative gene expression value. An unsupervised consensus clustering analysis was performed to identify edge perturbation-based clusters of GBM samples. Results revealed that the edge perturbation of GBM samples was stronger than that of normal samples. Four edge perturbation-based clusters of GBM samples were identified and showed prominent heterogeneities in prognosis, clinicopathological features, somatic genomic alterations, immune landscapes, and immunotherapeutic responses. In addition, a sample-specific perturbation of gene interaction score (SPGIScore) was constructed based on the differently expressed genes (DEGs) among four clusters, and exhibited a robust ability to predict prognosis. In conclusion, the bioinformatics approach based on sample-specific edge perturbation in gene interaction network provided a new perspective to understanding the multidimensional heterogeneities of GBM.

Establishment of a Necroptosis-Related Prognostic Signature to Reveal Immune Infiltration and Predict Drug Sensitivity in Hepatocellular Carcinoma.

Background: Hepatocellular carcinoma (HCC) is a common type of primary liver cancer and has a poor prognosis. In recent times, necroptosis has been reported to be involved in the progression of multiple cancers. However, the role of necroptosis in HCC prognosis remains elusive. Methods: The RNA-seq data and clinical information of HCC patients were downloaded from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. Differentially expressed genes (DEGs) and prognosis-related genes were explored, and the nonnegative matrix factorization (NMF) clustering algorithm was applied to divide HCC patients into different subtypes. Based on the prognosis-related DEGs, univariate Cox and LASSO Cox regression analyses were used to construct a necroptosis-related prognostic model. The relationship between the prognostic model and immune cell infiltration, tumor mutational burden (TMB), and drug response were explored. Results: In this study, 13 prognosis-related DEGs were confirmed from 18 DEGs and 24 prognostic-related genes. Based on the prognosis-related DEGs, patients in the TCGA cohort were clustered into three subtypes by the NMF algorithm, and patients in C3 had better survival. A necroptosis-related prognostic model was established according to LASSO analysis, and HCC patients in TCGA and ICGC were divided into high- and low-risk groups. Kaplan-Meier (K-M) survival analysis revealed that patients in the high-risk group had a shorter survival time compared to those in the low-risk group. Using univariate and multivariate Cox analyses, the prognostic model was identified as an independent prognostic factor and had better survival predictive ability in HCC patients compared with other clinical biomarkers. Furthermore, the results revealed that the high-risk patients had higher stromal, immune, and ESTIMATE scores; higher TP53 mutation rate; higher TMB; and lower tumor purities compared to those in the low-risk group. In addition, there were significant differences in predicting the drug response between the high- and low-risk groups. The protein and mRNA levels of these prognostic genes were upregulated in HCC tissues compared to normal liver tissues. Conclusion: We established a necroptosis-related prognostic signature that may provide guidance for individualized drug therapy in HCC patients; however, further experimentation is needed to validate our results.

Connexin Mutations and Hereditary Diseases.

Inherited diseases caused by connexin mutations are found in multiple organs and include hereditary deafness, congenital cataract, congenital heart diseases, hereditary skin diseases, and X-linked Charcot-Marie-Tooth disease (CMT1X). A large number of knockout and knock-in animal models have been used to study the pathology and pathogenesis of diseases of different organs. Because the structures of different connexins are highly homologous and the functions of gap junctions formed by these connexins are similar, connexin-related hereditary diseases may share the same pathogenic mechanism. Here, we analyze the similarities and differences of the pathology and pathogenesis in animal models and find that connexin mutations in gap junction genes expressed in the ear, eye, heart, skin, and peripheral nerves can affect cellular proliferation and differentiation of corresponding organs. Additionally, some dominant mutations (e.g., Cx43 p.Gly60Ser, Cx32 p.Arg75Trp, Cx32 p.Asn175Asp, and Cx32 p.Arg142Trp) are identified as gain-of-function variants in vivo, which may play a vital role in the onset of dominant inherited diseases. Specifically, patients with these dominant mutations receive no benefits from gene therapy. Finally, the complete loss of gap junctional function or altered channel function including permeability (ions, adenosine triphosphate (ATP), Inositol 1,4,5-trisphosphate (IP3), Ca2+, glucose, miRNA) and electric activity are also identified in vivo or in vitro.

An integrative multi-omics analysis based on liquid-liquid phase separation delineates distinct subtypes of lower-grade glioma and identifies a prognostic signature.

BACKGROUND: Emerging evidences have indicated that the aberrant liquid-liquid phase separation (LLPS) leads to the dysfunction of biomolecular condensates, thereby contributing to the tumorigenesis and progression. Nevertheless, it remains unclear whether or how the LLPS of specific molecules affects the prognosis and tumor immune microenvironment (TIME) of patients with lower-grade glioma (LGG). METHODS: We integrated the transcriptome information of 3585 LLPS-related genes to comprehensively evaluate the LLPS patterns of 423 patients with LGG in The Cancer Genome Atlas (TCGA) cohort. Then, we systematically demonstrated the differences among four LLPS subtypes based on multi-omics analyses. In addition, we constructed the LLPS-related prognostic risk score (LPRS) for individualized integrative assessment. RESULTS: Based on the expression profiles of 85 scaffolds, 355 regulators, and 3145 clients in LGG, we identified four LLPS subtypes, namely LS1, LS2, LS3 and LS4. We confirmed that there were significant differences in prognosis, clinicopathological features, cancer hallmarks, genomic alterations, TIME patterns and immunotherapeutic responses among four LLPS subtypes. In addition, a prognostic signature called LPRS was constructed for individualized integrative assessment. LPRS exhibited a robust predictive capacity for prognosis of LGG patients in multiple cohorts. Moreover, LPRS was found to be correlated with clinicopathological features, cancer hallmarks, genomic alterations and TIME patterns of LGG patients. The predictive power of LPRS in response to immune checkpoint inhibitor (ICI) therapy was also prominent. CONCLUSIONS: This study provided a novel classification of LGG patients based on LLPS. The constructed LPRS might facilitate individualized prognosis prediction and better immunotherapy options for LGG patients.

A Pyroptosis-Related Gene Prognostic Index Correlated with Survival and Immune Microenvironment in Glioma.

PURPOSE: As an inflammatory form of programmed cell death, pyroptosis has been well established to be associated with tumorigenesis and tumor immune microenvironment. In this paper, we aimed at the construction of a pyroptosis-related gene prognostic index (PRGPI) for predicting prognosis and guiding individualized immunotherapy in glioma patients. PATIENTS AND METHODS: Pyroptosis-related genes (PRGs) were identified based on a detailed review of published literatures. The transcriptome data and clinical information of glioma patients were obtained from CGGA and TCGA databases. PRGPI was constructed by using the multivariate Cox regression method. The immune cell infiltration level was analyzed via CIBERSORT algorithm. The tumor immune dysfunction and exclusion (TIDE) algorithm was applied to evaluate the potential response to immune checkpoint inhibitor (ICI) therapy. The expression patterns of PRGs in PRGPI were validated in cell lines and pathological specimens. RESULTS: We identified a total of 31 PRGs. Among them, PRGs (CASP3, DPP9, MAPK8, PELP1 and TOMM20) were selected for the construction of PRGPI. In both training (CGGA693) and validation (CGGA325 and TCGA) cohorts, PRGPI-high patients showed an inferior survival outcome compared with PRGPI-low patients. ROC curves illustrated that the prognostic prediction power of PRGPI was robust. A nomogram was developed based on independent prognostic indicators (PRGPI, age and WHO grade), and also exhibited a strong forecasting ability for overall survival (OS). Additionally, PRGPI-high patients exhibited higher immune, stroma and ESTIMATE scores, lower tumor purity, higher infiltration of M2-type macrophages, lower infiltration of CD8+ T cells and activated NK cells, higher tumor mutation burden (TMB), and higher expression of immune checkpoints. TIDE showed that PRGPI-high group had more responders of ICI therapy than PRGPI-low group. Finally, the expression patterns of five selected PRGs in PRGPI were significantly different between normal and glioma. CONCLUSION: The constructed PRGPI can be used for predicting prognosis and guiding individualized immunotherapy in glioma patients.

Identification of Critical m6A RNA Methylation Regulators with Prognostic Value in Lower-Grade Glioma.

Increasing evidences have revealed that N6-methyladenosine (m6A) RNA methylation regulators participate in the tumorigenesis and development of multiple tumors. So far, there has been little comprehension about the effects of m6A RNA methylation regulators on lower-grade gliomas (LGG). Here, we systematically investigated the expression profiles and prognostic significance of 36 m6A RNA methylation regulators in LGG patients from the TCGA and CGGA databases. Most of the m6A RNA methylation regulators are differentially expressed in LGG tissues as compared with normal brain tissues and glioblastoma (GBM) tissues. The consensus clustering for these m6A RNA methylation regulators identified three clusters. Patients in cluster 3 exhibited worse prognosis. In addition, we constructed an m6A-related prognostic signature, which exhibited excellent performance in prognostic stratification of LGG patients according to the results of the Kaplan-Meier curves, ROC curves, and univariate and multivariate Cox regression analyses. In addition, a significant correlation was observed between the m6A-related prognostic signature and the immune landscape of the LGG microenvironment. The high-risk group exhibited higher immune scores, stromal scores, and ESTIMATE scores but lower tumor purity and lower abundance of activated NK cells. Moreover, the expression level of immune checkpoints was positively correlated with the risk score. To conclude, the current research systematically demonstrated the prognostic roles of m6A RNA methylation regulators in LGG.

A Prognostic Ferroptosis-Related lncRNAs Signature Associated With Immune Landscape and Radiotherapy Response in Glioma.

Recent studies have demonstrated that long non-coding RNAs (lncRNAs) are implicated in the regulation of tumor cell ferroptosis. However, the prognostic value of ferroptosis-related lncRNAs has never been comprehensively explored in glioma. In this study, the transcriptomic data and clinical information of glioma patients were downloaded from TCGA, CGGA and Rembrandt databases. We identified 24 prognostic ferroptosis-related lncRNAs, 15 of which (SNAI3-AS1, GDNF-AS1, WDFY3-AS2, CPB2-AS1, WAC-AS1, SLC25A21-AS1, ARHGEF26-AS1, LINC00641, LINC00844, MIR155HG, MIR22HG, PVT1, SNHG18, PAXIP1-AS2, and SBF2-AS1) were used to construct a ferroptosis-related lncRNAs signature (FRLS) according to the least absolute shrinkage and selection operator (LASSO) regression. The validity of this FRLS was verified in training (TCGA) and validation (CGGA and Rembrandt) cohorts, respectively. The Kaplan-Meier curves revealed a significant distinction of overall survival (OS) between the high- and low-risk groups. The receiver operating characteristic (ROC) curves exhibited robust prognostic capacity of this FRLS. A nomogram with improved accuracy for predicting OS was established based on independent prognostic factors (FRLS, age, and WHO grade). Besides, patients in the high-risk group had higher immune, stroma, and ESTIMATE scores, lower tumor purity, higher infiltration of immunosuppressive cells, and higher expression of immune checkpoints. Patients in the low-risk group benefited significantly from radiotherapy, while no survival benefit of radiotherapy was observed for those in the high-risk group. In conclusion, we identified the prognostic ferroptosis-related lncRNAs in glioma, and constructed a prognostic signature which was associated with the immune landscape of glioma microenvironment and radiotherapy response.

Identification and Validation of a Prognostic Immune-Related Alternative Splicing Events Signature for Glioma.

BACKGROUND: Glioma is the most common malignant brain tumor in adults, with its tumor-promoting immune microenvironment always being intricate to handle with. Amounts of evidence has accumulated to suggest that alternative splicing (AS) is related to tumor immune microenvironment. However, comprehensive analysis of immune-related AS events and their clinical significance are still lacking in glioma. METHODS: AS events and transcriptome data of 653 glioma patients were downloaded online. ssGSEA was performed on transcriptome data of 653 patients to divided them into low, medium and high immune cell infiltration groups. Immune-related AS events were filtrated based on this grouping. Then lasso Cox regression analysis and multivariate Cox regression analysis were done to achieve an immune-related AS events prognostic signature for glioma. Kaplan-Meier analysis, ROC analyses, univariate Cox regression and multivariate Cox regression were performed to reveal the independent prognostic role of this signature. Meanwhile, a nomogram was constructed to achieved better prognostic value for glioma patients. Besides, functional enrichment analyses and correlation analyses with immune cells infiltration were used to validated the immune-related characteristic of this signature. RESULTS: 36 immune-related AS events were achieved based on the grouping mentioned above. A nine-immune-related alternative splicing event signature was built for glioma patients. This signature showed an independent prognostic value and a nomogram containing gender, age, Karnofsky performance score, grade, IDH status, MGMT promoter status and risk score derived from the signature was constructed with a higher predictive ability for overall survival. Association with the infiltration of immune cell subtypes was validated and functional enrichment analysis found that the signature was mainly enriched in immune-related and pro-tumor functions. CONCLUSION: Our research presented all immune-related AS events in glioma, identified an immune-related prognostic AS events risk model and a nomogram was constructed to predict the prognosis individually and more precisely. Tight connection was verified between this signature and clinical characteristics. Also, immune cells infiltration and immune checkpoints expression level were proved to link to risk scores, which enhanced the understanding of relationship between AS events and glioma immune microenvironment, firstly revealing the potential role of AS in immunotherapy of glioma.